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VIII. NONLINEAR OPTICS -- QUANTUM PICTURE: 45 

A Quantum Mechanical View of the Basics of Nonlinear Optics 46 

In what follows we draw on the discussion of the density operator in Review of Basic 
Quantum Mechanics: Dynamic Behavior of Quantum Systems, Section II of the lecture 
set entitled The Interaction of Radiation and Matter: Semiclassical Theory 
(hereafter referred to as IRM:ST). The macroscopic polarization is given by 

(i>)=Tr(pP) [VIII-1] 

We take the total Hamiltonian of a particular system in the form 

H = X + X mt +X rdndom [VIII-2] 

where 0-[ random is a Hamiltonian describing the random perturbations on the system by 
the thermal reservoir surrounding the system. Thus 



where 

dt 



dp 



[ VIII-3 ] 

relax 



= {[P'X -don,] [ Vffl-4 ] 

relax 



To find the nonlinear susceptibility, we make use of the following perturbation 
expansions: 

p= p (0) +p (l) +p (2) + ... [VIII-5a] 



See Nonlinear Optics — Classical Picture which is Section VII in the lecture set entitled On Classical 
Electromagnetic Fields (OCEF). 

See, for example, Chapter 2 in Y. R. Shen's Principles of Nonlinear Optics, Wiley (1984). 
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with 



/ p \_/ p(0) \, / p(D \ , / 5(2) \ 



(P (a) )=Tr(p (a) P) 



[ VIII-5b ] 
[ VIII-5c ] 



By substituting these expansions into Equation [ VIII-3 ] and equating terms of like 
order in H mX , we obtain the following hierarchy of equations: 



dp 



(o) 



Jo) 



'^°J + dt 



dp 



(1) 



dt 

d_p^ 

dt 



dt h 

i{[p<-U ] + [ P <»U w ]} + f 



relax 



= i{[p (2) ^o] + [p ( 'U,]} + ^ 



^ = -;{ [p( 3)^ o ] +[p(2Umt ]} + 



dp 



(3) 



relax 



relax 



relax 



[ VIII-6a ] 



[ VIII-6b ] 



[ VIII-6c ] 



[ VIII-6d ] 



Following earlier considerations, it is reasonable to write 



^dt 



l^=-r„„<«|p|«'>=-r„„-p„ 



relax 



[ VIII-7 ] 



Since the perturbing field is resolvable into an appropriate set of Fourier components 
(either discrete or continuous), we may write 
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# int = E# mt K) exp (-/ to,?) [ VIII-8a ] 

P=£p(o>*)exp(-i(D^) [VIII-8b] 

and resolve Equations [ VIII-6 ] into a hierarchy of algebraic equations. The first 
member of that hierarchy becomes 47 

pLW{^K--(oJ + r„„j=i<n|{p< )^ i J( J-^ i J % )p< )}|n') 

[ VIII-9a ] 



or 48 



piVK) = ^{«l^ i „,K)l«){p[?-piv}^K'-»,;r nn ') . [vm-9b] 



The second member of the hierarchy becomes 

P£'K){ ^'K'-«,)+r„„,} = ^i:<n|[p () (Q)^),^ mt (Q),-C0^) ]\n) [VHI-lOa] 



or 



" *' »" [VHI-lOb] 

>K>*-M|n>[VK)} 



we presume in this development that 9{ int is a Hermitian operator. 

Recall that ^ W — v; w) = [i ( W — v) + w] is the so called complex Lorentzian denominator. 
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Thus, we see that this expansion yields results akin to those embodied in Equations 
[VII-8 ] through [ VIII- 13 ] in OCEF. 

Two-PHoton Absorption 

Let us consider two-photon absorption by a single atom (or by a set of paired k -states 
in a semiconductor as discussed Section 7 of IRM:ST). 



C0 B 


vf0 2 


/CO! 





^\C0^ 




COb t 


C0 A 




p- 



Two-Photon Absorption Process 

Tfrom earlier discussions, the appropriate interaction Hamiltonian is given by 

^ 2f =-rr-E (+) (R,()=- e ?-[E! +) (R,f) + E; +) (R,| [VIII-11] 

In the appendix to this section we establish — viz., in Equation [ VIIIA-12 ] — that the 
second approximation to the transition rate is given by 

— =— V P {2) (tt) 

x (2) " dtr ' 



tr 



co,-co m 



5(co,. -(o f ) 



Thus for the two-photon process 
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r (2) 



2%e 4 



</|? E (+) |m)(m|r-E (+) |j) 



8(g),. -©,) [VHI-12] 



For two beam, two-photon processes 



rP) 



2 Tie 4 



{/[r.E( + )[m)<m[f.Ef)[/) 



+ 



{/Ir-E^lrnXmlr-E^IO 



(0 2 -C0, 



[VIII- 13a] 



8(a) 1 +a) 2 -a) / ) 



and for one beam, two-photon processes 



2 Tie 4 



r (2) 



I 



</|rE (+) |m)(m|r-E (+) |i) 



(0, -CD, 



5(2to beam -to / )[Vin-13b] 



If we use once again the factorization used in Section VI of this lecture set — viz. we 
write the initial state as | i ) = \i= {i'J"Y} = \ A (i) )| F{i") } where | A{i) ) and 
| F{i") ) are, respectively, the initial electronic (atomic) and field (photon) states. 

Using closure on the intermediate photon states, we find for two beam, two- 
photon processes 



xTr [pE'-'E'-' E 2 +) E< +) ] 



[VIII- 14a] 



where 
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M (1,2) s£ 



(A(/)|r|A ^J.e^j.^/n^rlAd-)) 



+ 



{A(/)|?|A(m))-e 2 e,-(A(m)|?|A(0) 



CO, -co 

2 m 



[VIII- 15a] 



and for one beam, two-photon processes 

1 2rce 4 



f , ,fW (1,1)| 8(2©-©,) 

Tr [pEj"' Ej"' E< +) E| +) ] 



[VIII- 14b] 



x 



where 



M (U) = E 

m 

Four Wave Mixing: 



<A(/)|r|A (/*))•£, e r (A(/n)|r|A(i)> 



CO,-CO m 



[VIII- 15b] 



Four Wave Mixing Spectroscopy: 

We explore here a particular set of four wave interactions. To that end, we consider an 
input field which consists of two plane waves 



E„(r,f)= E a (k 1 ,co 1 )exp[/(k 1 -f-co 1 f)] 

+ E a(k 2 < co 2 )exp[/(k 2 -?-co 2 f) 



[VIII- 16 ] 



+ c.c. 



where frequencies COj and C0 2 lie in the visible part of the spectrum. As Equation 
[ VIII- 3 ] of OCEF informs us, this input incident on a third-order nonlinear material, 
will directly generate a component of nonlinear polarization at a frequency 2c0[ - CO 2 

~ viz. 
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p; l, (M= 3%S; s K;g) i; -cq 2 ) Ep^.toJ E T k, ^ E 8 [k 2 . co 



x 



exp (2k,- k 2 ) ■ r| exp [-z (2co 1 -03 2 ) f] 



[VIII- 17 ] 



which may be associated with the destruction of two photons at (0 l and the creation of 
photons at co 2 and 2(0,-0)2. 

1 




Four wave mixing - direct process 

If, however, the material lacks a center of symmetry there is also an alternate path to 
obtain 2cfl I -Cfl 2 radiation. In a second-order material [ VIII- 3 ] of OCEF informs us 
we can also have a nonlinear polarization at 2(0, -(0 2 



(2) p: l, (M= {2xS>i-o> 2 ; mjEpfr-k* (o-tojE^to, 

+%5 ) Y (2a) 1 ;-a) 2 )E fi (2k 1 ,2a) 1 )E;(k 2 ,a) 2 )} [VIII-18] 
xexp[z (2kj- k 2 )-fj exp [-i (2(0,-0)2) t] 
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if we have driving fields Eplkj- k 2 , C0j-C0 2 J and E p [2k [; 2(0, J. Such fields are, in 
fact, generate by the input field through the nonlinear polarization components 

^Prfr-k^-co^x^ [VIII-19a] 



and 

(2) Pr(2 k, 2(0,) = 2 X < 2 3 > 1 ;to 1 ) Epfk.ODj E T (k r coj [ VIII-19b ] 




Four wave mixing - indirect processes 



We may solve for Epfkj - k 2 ,co, -C0 2 ) and E p [2k,, 2C0J from Equation [ VIII-23 ] in OCEF — viz. 

{|Q| 2 S aP - QaQp - K/c 2 e ) £ ap (to)} E p (Q, to) co 2 P a (NL) (Q, to) [ VIII-20a ] 
or using the operator §(k) = (k • k) 1 - k k = |k| 2 |l- k kj defined in OCEF 

{ e £" a ; (co)3 ps (Q) - (co 2 /c 2 )6 aS } E 5 (Q, to) = (co 2 /c 2 ) (co) P p (NL) (Q, co) [ VIII-20b ] 
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If we define 

^ a8 (Q;a))^{8 £- a ;(a))3 p8 (Q)-(a)7c 2 )8 a8 } [VDI-21 ] 

we can, in principle , write 

E a (Q, co)= [co 2 /c 2 ] § i (Q ; to) eg (to) P 8 (NL) (Q. to) [ VIII-22 ] 

Therefore, for a material lacking a center of symmetry, the complete nonlinear 
polarization may be written 

C L (f,r)= 3exp[/(2k 1 -k 2 )-f]exp[-/(2to 1 -(0 2 )?]Ep(k 1 ,to)E y (k 1 ,o) 1 )Ejk 2 ,to 2 ) 
x{XaPTB K;to i; -to 2 ) 

4 [VIII-23] 
+- [«7c 2 ] X^ ) Y (co 1 -co 2 ;co 1 )^J(k 1 -k 2 .o) 1 -o) 2 ) e^©,-©,) x$ (o) i; -(o 2 ) 



3 

2 



+ 



^ H^ 2 ] x^Kj-coJs -[(2k,, 2to,) 64(20)0 x^K;^)} 



Notice the "resonance" in the second term when |kj- k 2 | and oOj — to 2 are "tuned" to the wave number and frequency 
of coupled electromagnetic excitations or normal modes such as plasmons, polaritons and magnons. 



Degenerate Four Wave Mixing (DFWM) - Phase Conjugation: 

Consider now the reflected waves generated by the direct process — i.e., Equation 
[ VIII- 17 ] — when the input field is the sum of two counter propagating "pump" fields, 
with [k p ,to] and [-k p ,co] , plus a complex wave E,„(r, t) represented a range of wave 

numbers components |{k f },u)j --viz. 
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E„(f , t) = |E< +) exp [i (k p ■ r-(|> + )] + E<-» exp [i (-k p ■ f -<> _) 
+£e^ exp[/(k,, r-^^j lexp[-/co t] + c.c 



as shown below 



input 
signal 





/ 

/ nonlinear 
' crystal 



[ VIII-24 ] 



According to Equation [ VIII- 17 ], reflected waves are generated by nonlinear 



polarizations components with |{k f },3o)j, [{2 k p + k f },Q)j , ||2 k p - k f },-co 



, and 



. The latter component is the so-called phase conjugate term and may 



[{kJ,-G3 

be written 

(3) Pr(f,r)= ^jywcojE^ E< +) (k p , ©,) exp [-i«> + +<>_)] 



^Eg () exp[-/(k,, -r-^^j lexp[-/co t] + c.c. 



[ VIII-25 ] 



which may be interpreted as the time reversed image of E in (r, t) . 



TIME REVERSED IMAGING 
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i ncident signal 




The Raman Effect 

The Raman Nonlinearity: 

In our discussion of the dielectric susceptibility thus far we have been implicitly 
assuming that the nuclei of system were clamped into fixed positions and have focused 
on the nonlinearities associate with electron dynamics. This is too limited a view. In 
general, we may expect that polarizability of, say, a molecule is a function of the 
nuclear positions and that the dipole moment can expressed 



|u | denotes a collection of nuclear positions where u is the displacement of the s 

nucleus from its nominal equilibrium position. Further, we may expand these terms in 
powers of the displacement — viz. 



P„({«" ) },E)=P«"'({«" ) }) + a„ p ({u l '>}) 




[ VIII-26 ] 





(o) 



H) 



i 




a 



[ VIII-27a ] 
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■ J ({a (i) },E)^a ap ({a W })E p -agE„+£flW E p u« + - [ VIII- 27b ] 



The second term in Equation [ VIII-27b ] is the so called Raman nonlinearity. 

has the dimensions of charge, Equation [ VIII-27a ] is 



Since 



usually written 

Pl 0) ({fi ( ' ?) })-PL 0) +E4(^)< 5) + -- [VIII-27a'] 

where <?* p (.v) is the dynamic effective charge tensor associated with the s 

nucleus. The nuclear displacements are a linear combination of the normal modes of 
the system 

3N 

Up )= E r p (s;v)g(v) [VHI-28] 

v=l 

so that we may recast Equations [ VIII- 27a ] and [ VIII- 27b ] in the form 

PL 0) ({fi W }) = P?+I<(v) iv) + - [ VIII-29a] 



V=] 

N 3N 

P? = I>$y Ep nW =Eflff(v) E p q(v) [ VIII-29b ] 

i=l v=l 



AT iV 

where c*(v)= r p (s; v) and a® (v)=£a^ r y (s;v) 
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Quantum Theory of Raman Scattering 

The Raman nonlinearity provides a coupling between the photon field described by the 
Hamiltonian 



Xrad =E fc»S 



ka 



a r ar +- 

k,a k a ^ 



[ VIII-30 ] 



and the phonon field described by the Hamiltonian 



v r + i 



[ VIII-31 ] 



where 6 v and £ v are, respectively, the phonon creation and destruction operators. In 
the electric dipole approximation, the interaction Hamiltonian can be written 



-. 3N 

^,„ t =-2W( v ) E « E P <?( v ) 



v=l 



[ VIII-32 ] 



In earlier discussions we have shown that 



q (v)- 



f * V/ 2 
n 



2M(0, 



[ VIII-33a ] 



and 



kc' 



1/2 



V 



[VEI-33b] 



so that the complete interaction Hamiltonian may expressed as 
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2% h 



3;Y 



h 



.1/2 



v=l V 



2Mm 



[ VIII-34 ] 



ka 



k', o' 



To concentrate on the key issue, we write scattering Hamiltonian which is 
proportional to a and which describes the scattering of a photon from the 

|k',o'} (initial) state to the {k,G } (final) state 



where a {R) (l,F; v) = aff(v) e a (/) e p (F) 



271 ft 



.1/2 



2M0) 



(to, C0 F )' /2 a (R) (/,F; v) a J a, [6 1 + 6 v ) [ VIII-35 ] 



v y 



We may now use the Fermi golden rule 49 to calculate the transition rate per 
"molecule" for inelastic Stokes scattering events -- viz. for 
{n n n F ,n v } -^{n I -\,n F + \,n v + \} 



2k 

2 



n 



2% h 
v 



2Mm 



(to, ® F ) |a (R) (/^;v)f 



[ VIII-36a ] 



x|(/i / — + + dpdj &l \nj,n F ,n v )\ 8((0 / -Q) F -C0 v ) 



1 8n- 



r(S) 



V 7 



2Mco 



v ; 



(t0 / C0 F )| a [R) (l,F; v)f n, (n F + l) (n v + l) 8(a), -0) F -G) v ) [ VIII-36b ] 



See the appendix to this section for a recapitulation of the theory of transition rates. 
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and the transition rate per "molecule" for inelastic anti-Stokes scattering events 
viz. for {n,, n F , t%] -l,n F + l,r\, -1} 



1 2% 



x< AS > h 2 



8tt 



2% h 



n 



2M 0) 



V J 



(to, co F )|a (R) (/,F;v)| 2 



y.\n I -\,n F + \,n v -\\ a\ a, b v |n, , n F ,n v f 8(a), +C0 v -C0 F ) [ VIII-36c ] 



V 1 



2M(0 



v 7 



(to, co F )|a (R) (/,F;v)| 2 ra, (n F +\) n v 8(a), +to v -co F ) 



For a total number of "molecules" Vp mol we may write the total time rate of change of 
photons of frequency G) F as 



d n. 



d t 



d n F 



8?1 Pjnoi t«/ 1«/ 

y 



V 2M co v , 



a (R) (/,F;v) 



Xn ( (n F + l) (n v + l)8((Dj -(0 F -(0 V ) 



d t 



A S 



y 



V 2M oo v , 



« (R, (/^;v) 



xn, (n F + l) n v 8((0 7 +CO v -CO F ) 



[ VIII-37a ] 



[ VIII-37b ] 



Stimulated Raman Effect: Simple model 

Suppose that a high power laser at (Q, — the pump -- propagates parallel to the z-axis 
and that the Raman medium occupies the half-space z > 0. The spatial rate of change 
of the Stokes photons at C0 F — the signal -- may be written 



50 



^=<Wz)Mz)+i] 

dz 



[ VIII-38 ] 



Assuming for the nonce that n y ~ . 
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where 



8rc 3 P mo r©/ CO, 



cV y. 



n 



2M co„ 



a (R) (l,F;v) 



[ VIII-39 ] 



v V^ J " ™v J 



0.20 



40 W/crrrt2 
60 W/cm*2 
100 Wtei7f^2 
300 Wfcm*2 




-10 -5 5 10 15 20 25 30 

Probe Detuning (MHz) 

Raman Gain in rubidium vapor for different coupling laser intensities 



The rate constant y v is included this expression for gain to account for the broadening 
of a given vibrational state. Each Stokes event, of course, also depletes the incident 
beam so that 



a z 



[ VIII-40 ] 
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which means that n F (z) + n l (z) = n 1 (0) . [ VIII-41 ] 

Substituting this conservation condition into Equation [ VIII-38 ] we see that 

^M =Gr [ W/ ( )_ n F (z)] [n F {z) + l] [ Vm-42a ] 

clz 



or that 



Therefore 



d~ nJz) dnJz) r / o 

r - ( - r ( \ -m = G « 1+ n > d z [ VIII " 42b ] 



exp{ G R [1 + ^(0)] z}-l 

niz) = n 7 - r . . [ VIII-43 ] 

' A \(0) + exp{G R [l + n 7 (0)]z} 

This simple model predicts that Stokes photons initially increase linearly with distance 
and that there complete conversion at distances large compared to j G R [l + n 7 (0)] j 



i.e. 



lim n F (z) = G R [l + n, (0)] z [ VIII-44a ] 

]imn F (z) = n I (0) [ VIII-44b ] 
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n,(0) = l0 




0.4 0.8 1.2 1.6 



Coherent anti-Stokes Raman Spectroscopy (CARS): 

There is an important four wave mixing process involving the Raman nonlinearity. We 
may see how this comes about by realizing that Equation [ VIII-29 ] implies a potential 
energy contribution 

«£flS ) (v)E B E | ,9(v) [VIII-45 

1 v=l 

so that inclusion of the Raman nonlinearity modifies the equation of motion for the 
vibrational normal modes - viz. 

q(t,v)+0>l q(t,v)= ± «5'(V) E a E p [ VIII-46 ] 

When the normal mode is driven by a two (visible) wave input (see discussion of four- 
wave processes above), its dominant response will be at the difference frequency 
03j-G) 2 — viz. 
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q(t,v) +(0 2 v q(t,v) 

1 [ VIII-47 ] 

= — «S'( V ) E >i) Ep(t0 2 )exp[-j {(Q-(Q 2 ) t\ + c.c. 



so that 



9 (?,v) = — p 7 , — ; 2 — — [vm-48] 

From Equation [ VIII-29b ] we see that excitation of a vibrational mode, through the 
Raman nonlinearity, generates a nonlinear polarization 



1 v a V( v ) «i ( v ) E pN E tN e sK) exp -i (20)-© 2 ) f] 
P a = — 2- — ; n [ vm-49 ] 

M v=1 to'-^-COj 



3iV 



at the four wave mixing frequency 2c0j-C0 2 . Thus, we see that the medium has an 
effective third order susceptibility 

[Zg,(2^-»J] =X^+^4t^ [ Vm - 5 °] 
which exhibits a resonance when CO, - CO 2 equals C0 V . 
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